Fluorescence decay in aperiodic Frenkel lattices 
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We study motion and capture of excitons in self-similar linear systems in which interstitial traps 
are arranged according to an aperiodic sequence, focusing our attention on Fibonacci and Thuc- 
Morse systems as canonical examples. The decay of the fluorescence intensity following a broadband 
pulse excitation is evaluated by solving the microscopic equations of motion of the Frenkel exciton 
problem. We find that the average decay is exponential and depends only on the concentration 
of traps and the trapping rate. In addition, we observe small-amplitude oscillations coming from 
the coupling between the low-lying mode and a few high-lying modes through the topology of the 
lattice. These oscillations are characteristic of each particular arrangement of traps and they are 
directly related to the Fourier transform of the underlying lattice. Our predictions can be then used 
to determine experimentally the ordering of traps. 

PACS numbers: 71.35+z; 36.20.Kd; 78.90.+t 



I. INTRODUCTION 

The discovery of quasicrystals |jj and the fabrication 
of aperiodic low-dimensional systems have encouraged 
the study of the physical properties of solids with long- 
range order, lacking translational symmetry. Physical 
aperiodic systems, whose structural order is described 
by means of deterministic sequences, present rather ex- 
otic electronic properties not shared by crystalline and 
amorphous solids. Peculiar signatures of these unique 
electronic states are Cantor-like spectra and self-similar 
wave functions. Several electronic properties of solids can 
be inferred from optical characterization techniques like 
optical absorption, photolumincscence and fluorescence 
after pulse excitation. Therefore, a complete understand- 
ing of the interplay between the electronic properties and 
the underlying aperiodic order requires a detailed anal- 
ysis of the optical dynamics of this kind of ordering of 
matter. This interest has motivated various works deal- 
ing with optical properties of aperiodic systems, mainly 
devoted to Fibonacci and Thue-Morse semiconductor su- 
perlattices. 

Following a long term project regarding electronic and 
transport properties of aperiodic systems, j?]-|lO| our 
group focused its attention on optical absorption spec- 
tra of Frenkel excitons in Fibonacci and Thue-Morse lat- 
tices. [ p]Jl^1 Our main aim was to learn about the phe- 
nomenology of incoherent and coherent exciton dynamics 
in molecular aggregates and polymers exhibiting long- 
range correlations. We considered Fibonacci and Thue- 
Morse lattices as canonical aperiodic models that are nei- 
ther periodic nor random. In turn, they are inherently 
close to realistic aperiodic systems where the effects of 
long-range order might change dramatically the exciton 
dynamics. We found numerically several characteristic 
lines specific of each aperiodic system that are not present 
in the spectra of periodic or random Frenkel lattices, 
thus being an adequate way for determining the particu- 
lar structural order of the system from experiments. pd| ] 
Moreover, using a model based on the Pauli master equa- 



tion, we have also studied motion and capture of incoher- 
ent excitons when traps are arranged according to aperi- 
odic sequences. p2j As an interesting feature, amenable 
to experimental confirmation by means of luminescence 
studies at moderate temperature, we obtained that the 
decay of the survival fraction of incoherent excitons is 
simply exponential P(t) ~ exp(— At) instead of the 
asymptotic stretched exponential P(t) ~ exp(— At 1 / 2 ) 
appearing in one-dimensional random lattices. |l3| ] It is 
well-known that the exponent of the stretched exponen- 
tial varies when passing from trapping of classical (in- 
coherent) to quantum (coherent) excitons in random lat- 
tices. [[14| Then a natural question arises in the context of 
aperiodic systems, namely whether the simple exponen- 
tial behavior of the survival fraction holds in the quantum 
case. 

In this work we investigate the decay fluorescence due 
to Frenkel excitons in linear systems with interstitial 
traps arranged according to the Fibonacci and Thue- 
Morse sequences. We have focused on the analysis of ape- 
riodic arrangements of traps in systems where optically 
active centers are placed in a straight line. Nevertheless, 
the model can also be applied to other possible geometric 
configurations of interest as, for example, zig-zag chains 
whose sp 2 hibridization is explicitly considered. To this 
end, we make use of a general treatment which allows 
us to study the dynamics of Frenkel excitons in these 
lattices, solve the microscopic equations of motion and 
find the fluorescence intensity as a function of time. The 
main results of this paper are the following. First, we 
find out the decay law of the fluorescence intensity at 
low temperature, and compare it to the above mentioned 
simple exponential trend at moderate temperature. Sec- 
ond, and most important from an experimental point of 
view, we obtain distinctive features of the decay curves 
P(t) allowing us to unveil the ordering of traps present 
in the system. By means of an analytical approach we 
are able to explain the origin of these characteristics os- 
cillations, which are related to the Fourier transform of 
the arrangement of traps. Thus we successfully relate an 
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optical property (fluorescence decay) with an structural 
one (topology of the lattice). 

The remaining of the paper is organized as follows. In 
Sec. II we describe our model Hamiltonian and the way 
we arrange the traps following the Fibonacci and Thuc- 
Morse sequences. Afterwards, in Sec. Ill we focus our 
attention on the fluorescence intensity following a broad- 
band pulse excitation and solve the microscopic equation 
of motion. We find numerically that, superimposed to 
an overall exponential decay, there exist small-amplitude 
oscillations characteristic of each particular system. In 
Sec. IV we develop a homogeneous medium approxima- 
tion which accounts for the exponential decay, and, in 
Sect. V we introduce a perturbative treatment that suc- 
cessfully explains the small-amplitude oscillations. In 
this way, we clearly demonstrate that the frequency of 
the oscillations depends on the specific aperiodic order of 
the distribution of traps and, in addition, that these os- 
cillations arise as a consequence of the coupling between 
the low-lying mode and few high-lying modes through 
the topology of the lattice. Section VI concludes the pa- 
per with a brief summary of results and some general 
remarks on their physical implications. 

II. PHYSICAL MODEL AND THEORY 

A. Model Hamiltonian 

We consider a system of N optically active two-level 
centers, occupying positions r„ on a linear regular lat- 
tice with spacing a. For our present purposes we neglect 
all thermal degrees of freedom. Therefore, the effective 
Frenkcl Hamiltonian describing this system can be writ- 
ten (we use units such that fi = 1) 

Here and a n are Bose operators creating and anni- 
hilating an electronic excitation of energy V n at site n, 
respectively. J n i (n ^ I) is the intersite interaction of 
dipole origin between centers n and I. In the case where 
all centers have equal and parallel dipole moment the in- 
teraction can be cast in the form J„/ = — J(a/\r n — r;|) 3 , 
where J is the coupling between nearest-neighbor cen- 
ters. In this paper J is the unit of energy and J -1 the 
unit of time. Moreover, since J n i is a rapidly decreasing 
function of the distance between centers, we will omit 
interactions beyond nearest-neighbors. 

According to previous works, |p^|-[l8|] the non- 
Hermitian term — iy~] n r n ajja n accounts for the irre- 
versible trapping of the exciton due to interstitial traps 
and the subsequent photon emission. Here r„ = T if the 
n center has a trap associated with it and r„ = other- 
wise, T being the trapping rate. The site energy V n can 
be written, without loss of generality, as (V) + AV n . The 



first term stands for the average excitation energy and 
we drop it by setting an appropriate reference. The sec- 
ond term is a random variable reflecting the offset energy 
due to the influence of the surrounding medium. Since we 
are mainly interested in the trapping dynamics, we ne- 
glect the effects of diagonal disorder (AV n — 0), noting 
that they could be easily incorporated in our treatment 
if necessary. 

B. Fibonacci lattice 

In this work we will be concerned with the Frenkel 
exciton dynamics and trapping processes in systems pre- 
senting an aperiodic distribution of interstitial traps. The 
Fibonacci lattice (FL) is the archetypal example of de- 
terministic and quasiperiodically ordered structure. Any 
arbitrary Fibonacci system presents two kind of building 
blocks. In our case, we choose those blocks as individ- 
ual centers with and without a trap associated with it, 
hereafter called T and S, respectively. The Fibonacci 
arrangement can be generated by the substitution rule 
S — > ST, T — > S. In this way, finite and self-similar 
aperiodic lattices are obtained by n successive applica- 
tions of the substitution rule. The nth generation lat- 
tice will have N = F n elements, where F n denotes a 
Fibonacci number. Such numbers are generated from 
the recurrence law F n = F n -\ + F n —2 starting with 
Fq = F\ = 1; as n increases the ratio F n _i/F n converges 
toward r = (v5 — l)/2 = 0.618 . . ., an irrational num- 
ber which is known as the inverse golden mean. There- 
fore, lattice sites are arranged according to the sequence 
STSSTSTS... This type of distribution of traps will 
be referred to as minority lattice since the fraction of T- 
centers (1 — t for large enough FLs) is smaller than the 
fraction of S'-center (r for large enough FLs). Moreover, 
it is worth noticing that T-centers appear isolated in the 
minority lattice. In the same way, one can generate the 
majority lattice by replacing S <-» T, thus obtaining the 
sequence T STT ST ST . . . In this case T-centers can 
appear isolated or paired. 

C. Thue- Morse lattice 

The Thue-Morse lattice (TML) is also constructed 
from two building blocks S and T. The substitution rule 
in this case is S — > ST, T — > TS and the nth genera- 
tion lattice contains 2" centers. It is clear that the frac- 
tion of T-centers is exactly 1/2 for any generation lattice. 
The resulting sequence STT ST S ST . . . is also aperi- 
odic and self-similar but not quasiperiodic. Notice that 
both S- and T-centers appear isolated or paired. The 
TML present an inversion center so that the lattice is 
invariant under the transformation S <-> T. 
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D. Computation of physical magnitudes 

Having presented our model we now briefly describe 
the method we have used to calculate the fluorescence 
decay following a short-pulse excitation. The calcula- 
tion involves the total dipole moment operator T> = 
Yin ^«( a n + a «): where d n is the dipole moment of the 
n center. Notice that we are restricting ourselves to the 
case of systems whose length is much smaller than the 
optical wavelength. After the excitation, the state vector 
at time t evolves according to \ip(t)} = exp(— iHt)T>\vac) , 
| vac) being the exciton vacuum state. Besides radiative 
effects, the fluorescence intensity, normalized to the value 
at t = 0, is given by [M 



P{t) = 



IWQ)I^))I 2 



(2) 



The fluorescence intensity can also be viewed as the prob- 
ability of finding an exciton in the k = mode at time 
t. 

A reliable method for determining P(t) numerically 
has been given by Huber and Ching. [pl| These authors 
introduced a set of correlation functions 



G n (t) = Y, dl 
i 



ya.c\a n (t)a\ | vac) 



(3) 



where a n (t) — exp(iH.t)a n exp(— iTit) is the annihilation 
operator in the Heisenberg representation. The function 
G n (t) obeys the equation of motion 



.dG n {t) 
dt 



H n iGi(t), 



(4) 



with the initial condition G„(0) = d n . The diagonal el- 
ements of the tridiagonal matrix H n i are — iT n whereas 
off-diagonal elements are simply given by — J. The micro- 
scopic equation of motion is a discrete Schrodinger-like 
equation on a lattice and standard numerical techniques 
may be applied to obtain the solution. Once these equa- 
tions of motion are solved, the fluorescence intensity is 
evaluated from the relationship 



P(t) = 



ISn d* n G n {t)Y 

IE„ KI 2 I 2 



(5) 



III. NUMERICAL RESULTS 

We have solved numerically the equation of motion (ji|) 
using an implicit (Crank-Nicholson) integration scheme. 
|l9f We note that energy is measured in units of J and 
time in units of J -1 and, since 4 J is the exciton band- 
width, the energy and time scales can be deduced from 
each particular experimental situation. Aperiodic lat- 
tices are generated using the inflation rules discussed 



above. In order to minimize end effects, spatial pe- 
riodic boundary conditions are introduced in all cases. 
The maximum integration time and the integration time 
step were 20 and 2 x 10~ 3 , respectively. Smaller time 
steps give essentially the same results. In most cases 
of practical interest, one finds that the trapping rate is 
much less than the exciton bandwidth, that is, T <C 4J. 
Thus, as typical values, we have studied lattices with 
r/J = 0.01, . . . , 0.05. The maximum lattice size un- 
der consideration was TV = F\2 = 233 for FLs and 
N = 2 s = 256 for TMLs. Finally, local dipole moments 
are taken to be d n = 1/ vTV. 
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FIG. 1. Tiritt? decay of fluorescence intensity 1 in (a) Fi- 
bonacci lattices with TV ^ F12 — 233 clnters^flid a frac- 
tion of traps c = 0.618 . . . (majority lattice, solid lines) 
and c = 0.382... (minority lattice, dashed lines), and (b) 
Thue-Morse lattices with TV = 2 8 = 256 centers. In all cases 
the trapping rate is, from top to bottom, V = 0.01, . . . , 0.05 
in units of J. 

In our numerical simulations we have found that the 
fluorescence intensity decays exponentially, for both FLs 
and TMLs, in the full time interval considered by us. 
Figure |l| shows the results for majority and minority FLs 
with TV = 233 and TMLs with TV = 256, but we have 
checked that this behavior is independent of the system 
size. As expected, the larger the trapping rate, the faster 
the decay of the fluorescence intensity. Moreover, for the 
same value of the trapping rate, this decay is faster on 
increasing the fraction of traps (minority FLs, TMLs and 
majority FLs). Hence, as occurs in the case of completely 
incoherent excitons, the decay of Frenkel excitons in 
FLs and TMLs resembles that of periodic lattices instead 
of random ones, where the decay is stretched exponen- 
tial. The average decay of the fluorescence intensity fits 
very well to an exponential of the form exp(— 2cFt) in 
the considered time interval, c being the fraction of traps 
in each lattice. This time dependence has also been con- 
firmed by fitting the data to stretched exponentials of the 
form exp(— At@), obtaining in all cases j3 — 1.0000±10~ 4 . 
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A detailed analysis of data plotted in Fig. |] shows the 
occurrence of small-amplitude oscillations superimposed 
to the average exponential decay of the fluorescence in- 
tensity The functional form P(t) in all cases considered 
can be expressed as 



7Yt 



P(t) = exp(-2cr<) [1 + T 2 p(t)] 



(6) 



where p(t) is an oscillatory function which is inde- 
pendent of the trapping rate. This function can 
be readily extracted from numerical data as p(t) ~ 
T~ 2 ln[P(t) exp(2cFt)]. Figure ^ shows the typical form 
of the p(t) function for anyone of the trapping rates T 
considered in our study. Remarkably this function is, on 
the one side, independent of the number of active centers 
present in the lattice and, most interestingly, it does not 
depends on the concentration of traps, i.e. it is exactly 
the same for both majority and minority versions of FLs 
— this invariance is, of course, trivial for TMLs — . Fi- 
nally, by comparing Figs. @(a) and (b) we realize that 
the oscillatory pattern of pjt) is characteristic of the ar- 
rangement of traps, in the sense that its general features 
are very sensitive to the kind of aperiodic order being 
considered. 
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FIG. 2. Plot'tffp^t) £ l'^ 2 lri[^m exp(2bl't)] as 1 a function 
of time for (a) Fibonacci Mid (b) Thue-tilxjrse lattice. Solid 
and dashed lines indicate numerical results and perturative 
prediction, respectively. Results are independent of the num- 
ber of lattice sites and the concentration of traps. 
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(r„/r - c)a* n a r , 



(7b) 



As a first approximation, we neglect the topology of the 
lattice and take H. ~ Hhm- Physically this implies re- 
placing the aperiodic distribution of traps by a uniform 
distribution with a renormalized trapping rate cr. To 
proceed we notice that the overlap between the state vec- 
tor at time t and the initial state vector may be calculated 
from Poll 



(V(0)|V(*)> = J2 l(fc|2?|vac)| 2 expHi^), 



(8) 



where \k) denotes the one-exciton eigenvectors of ( |7a| ) 
and Ek the corresponding eigenvalues. The Hamiltonian 
Hhm with periodic boundary conditions can be exactly 
diagonalized, yielding the eigenvectors [M 



ifc> = (£) ' E ex p( 2 ™^) a « 



vac) 



(9a) 



with k — 0, . . . , N — 1. The corresponding eigenvalues 
read 



E k 



-icT - 2 J cos 2tt 



N 



-icT + cjfc. 



(9b) 



Therefore, the dipole moment matrix elements in the 
main field approach satisfy the relation (fc|X>|vac) = 
Sko- Consequently one obtains that (ip{0)\tj;(t) ) = 
exp(— cTt — iu>kt) so that we obtain the exponential be- 
havior Phm(^) = exp(— 2cTt). This result shows that, 
within the tight-binding approximation, the overall ex- 
ponential decay of the fluorescence is a general property 
of ordered systems described in terms of the Hamilto- 
nian (Q), with independence of this order being periodic 
or aperiodic. We have confirmed this point by consider- 
ing periodic arrangements of traps with different basis. 
In all cases we have observed that the exponential be- 
havior predicted by the homogeneous medium approxi- 
mation treatment successfully describes the overall decay 
obtained by numerical simulations. On the other hand, 
the exact result obtained within the that approximation 
provides us with a ground starting point which will allow 
us to explain in full detail the small-amplitude oscilla- 
tions shown in Fig. ^. 



IV. HOMOGENEOUS MEDIUM 
APPROXIMATION 



V. PERTURB ATIVE APPROXIMATION 



In this section we develop a theoretical approach to 
explain the overall exponential decay of the fluorescence 
in aperiodic lattices. To this end we rewrite the system 
Hamiltonian TL = 7Yhm — iTH.T , where 

Hkm = -icT ^ a i a n - J ^ {al a n+i + a n+i a n), (7a) 



The replacement of a given aperiodic distribution of 
traps by a uniform one with a renormalized trapping rate 
is a crude representation of the real system. To develop 
a better approximation one must explicitly consider the 
effects of the topology on the exciton decay through the 
term — iVH^. To carry out such an approximation we 
have inspired in the nondegenerate perturbation theory 
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developed by Huber to study random systems with sub- 
stitutional traps. [^O) Let \k) the perturbed eigenvector 
describing the exciton state when the influence of the 
term 7ix is taken into account. To evaluate the per- 
turbed dipole moment matrix element, we expand \k) in 
the basis of the unperturbed eigenvectors \k) given in 



\k) = \k) 



^ E k - E, 1 



l^k 



whereas the perturbed eigenvalues are given by 
E k =E k -iT{k\H T \k), 



(10a) 



(10b) 



with E k being the unperturbed eigenvalues. Making use 
of the fact that the dipole moment matrix element in- 
volving unperturbed eigenvectors is (fc|2?|vac) = Sko, the 
overlap between \ip(Q)) and \ip(t)} can be readily found, 
leading to 



P(t)=PHM{t) 



i + r 2 z(t) 



1 + T 2 z(0) 

= Pau(t) {1 + 2r 2 Re[z(i) - z(0)]} + O (r 4 ) , (11) 



with 



\(0\Hi 



■exp(-iuj ko t), 



(12) 



All 



where we have defined ujko = u> k — ujq — 4J sm 2 (irk/N). 
To obtain the above result we have made use of the re- 
sult (k\H.T\k) = 0, which according to (10b), indicates 



that E k = E k in the first order perturbative treatment. 
Hence, recalling Eq. (p), one gets 



P(t) 



\(Q\H T \k)f 



(coswfcoi - 1) 



(13) 



All 



The most remarkable fact of this expression is that it 
relates the small-amplitude oscillations of the fluores- 
cence decay, an experimentally measurable magnitude, 
with the Fourier transform of the distribution of traps 
describing the topological ordering of the system. To 
demonstrate this point, we recall that unperturbed eigen- 
vectors are orthogonal. Thus 

(0|H T |fc) = (0|]T ^ata„|fc> 

n 

4^r xp K)- (14) 

which is nothing but the afore mentioned Fourier trans- 
form. A comparison of the perturbative prediction ( |l3| ) 
with the numerical result is given in Fig. where an 
excellent agreement is achieved for both Fibonacci and 
Thue-Morse aperiodic arrangements. At this point we 



should stress that the above result is valid for any arbi- 
trary lattice. It relates the Fourier pattern of the distri- 
bution of traps with the finer details of the fluorescence 
decay. Since each particular arrangement of traps leads 
to a different Fourier pattern of p(t), we arrive at the con- 
clusion that the small-amplitude oscillations characterize 
every distribution of traps in much the same way as the 
usual X-ray diffraction analysis characterize the distribu- 
tion of scattering centers in the space. It is also worth 
mentioning that it is not difficult to show from ( |l4|) that 
dual lattices, namely those lattices obtained from the re- 
placement T <-> S, present the same value of | (0 1 7Yt | fc) | 2 - 
In particular, minority and majority Fibonacci lattices 
should display exactly the same oscillations, as it is in- 
deed the case in our numerical results. 

In order to get a deeper insight into the small- 
amplitude oscillations in the case of aperiodic arrange- 
ment of traps, we relate our results with those reported 
in the various works dealing with the Fourier transform of 
these sequences, most of them in connection with X-ray 
diffraction and Raman scattering in aperiodic semicon- 
ductor superlattices, and neutron scattering in quasicrys- 
tals. p2|~pr| Let us consider, in the first place, the case 
of the FL, which is known to exhibit a pure point Fourier 
spectrum displaying well defined Bragg peaks. |2j| The 
Fourier intensity |(0| 7^x1^)1 2 of a FLs consists of a series 
of peaks located at values of the momentum 2Tr(k/N) of 
the form 27r(m + na), where m and n are two arbitrary 
integers and a = (\/5 + l)/2 is the golden mean. Each 
peak of the Fourier intensity leads to a large contribu- 
tion to the sum ( |l3| ) and, consequently, to a well-defined 
frequency in the oscillatory pattern of p(t). A simple in- 
spection of Fig. ||(a) indicates that the pattern of p(t) is 
dominated by a single frequency and, in fact, the Fourier 
analysis of the data gives the approximate value ~ 3.456 
for this main frequency. The frequency of the particular 
mode k p that significantly contributes to the oscillatory 
pattern is then given by ujk p — loq — 3.456, leading to 
k /N ~ 0.3798, the ratio being independent of the sys- 
tem size. Therefore the mode is located close to the top 
of the excitonic band (k/N = 0.5). 

Figure ||( a) shows |(0|HT|fc)| 2 /w 2 as a function of the 
ratio k/N evaluated from expression (|l4]), the results be- 
ing independent of the system size. A sequence of Bragg 
peaks, arranged according to successive powers of r, are 
clearly observed in this plot. The stronger one, which 
is responsible of the highest-frequency oscillations, is lo- 
cated at k/N ~ 0.38 ~ r 2 . Making use of the relation- 



ships T 



1 and a — 1 + r, we find the two indices 



labeling this peak to be m = 2 and n = — 1. The Fourier 
pattern is symmetric around the top of the excitonic band 
and, in fact, there exists also a similar peak located at 
k/N ~ 0.62 ~ 1 — t 2 = r, whose indices are m = — 1 
and n = 1. The remaining, in general weaker, peaks also 
give rise to oscillations in p(t) but their contribution is 
not appreciable in the decay curve because of they are 
closer to the bottom of the excitonic band and then lead 
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to low-frequency oscillations. 
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FIG. 3. (a) Iftot bf Ku^T^I^^'as a functioh of the ra- 
tio k/N for FLs sffowing trie occurrence of weU-defmed Bragg 
peaks, labeled according to succfes'rve powers of the inverse 
golden mean r. The arrow indicates the peak responsible of 
the largest-frequency oscillations observed in Fig. 2. (b) Plot 
of |(0|HT|fc)] 2 / w fco as a function of the ratio k/N for a TML 
with 2 8 = 256 active centers. A dense arrangement of peaks 
labeled according to the formula (2n + l)/(3 x 2 m ) is clearly 
observed. 

Let us consider now the case of the TML, which is 
known to exhibit a singular continuous Fourier trans- 
form, lacking well defined Bragg peaks. J2(| For finite 
TMLs the Fourier intensity consists of a dense arrange- 
ment of peaks whose number depends on the sample size. 
For a TML composed of 2" building blocks of the same 
size one expects to find 2" _1 peaks in the related Fourier 
transform whose positions are located at values of the 
momentum 2-K(k/N) of the form (2n + 1)/(3 x 2 m ), where 
n and m are two arbitrary integers. |26| ] An inspection of 
the Fig. |2|(b) indicates that, conversely to the case of FL, 
the pattern of the oscillatory fluorescence decay corre- 
sponding to the TML cannot be described in terms of 
a single main frequency. In fact, the Fourier analysis of 
numerical data shows a considerable number of frequen- 
cies closely packed. Figure ||(b) shows | (0|Ht|&) | 2 / w fco 
as a function of the ratio k/N for a TML with N = 2 a 
active centers. A dense distribution of peaks is clearly ob- 
served. We have checked that most peaks can be labeled 
by means of the formula (2n + 1)/(3 x 2 m ) . This we show 
in the figure for the most relevant peaks. By comparing 
Fig. ^(b) with X-ray results obtained for semiconducting 
TM superlattices, [^6| we see that the study of the oscil- 
lations in the fluorescence decay allow us to analyze the 
topological structure for a wide variety of aperiodic sys- 
tems, even if those system do not give rise to well defined 
Bragg peaks in the reciprocal space. 

Our results are suitable for a direct physical interpre- 



tation. Small-amplitude oscillations of largest frequency 
LUk p — loq are caused by the coupling of two modes, namely 
the lowest-lying and k p , through the topology of the ape- 
riodic distribution of traps. Notice that different arrange- 
ment of traps would lead to completely different Fourier 
intensity |(0|7iT|fc)| 2 and p(t), as it can be seen from a 
comparison of Figs. ||(a) and (b). Thus, the exciton acts 
as a probe of the spatial distribution of traps. This is one 
of the main results of the present work since it provides 
us with a useful method to be applied in experimental 
situations. To elucidate the particular arrangement of 
traps in a linear lattice, first we can analyze the fluo- 
rescence decay after (5-pulse excitation. If this decay is 
stretched exponential, the arrangement of traps is com- 
pletely disordered, whereas simple exponential behavior 
indicates absence of disorder corresponding to either pe- 
riodic or aperiodic distribution of trapping centers. The 
slope of the average exponential decay permits to find out 
the concentration of traps through the exponent value of 
PaM(t), and then p(t) ~ ln[P(t) / Phm^)]- By performing 
the Fourier analysis of p(t) one can determine its related 
power spectrum and then, comparing the obtained re- 
sults with theoretical predictions from ( |l3| ) and (|l4|), to 
extract relevant information about the topological order 
of the underlying structure. 

On the other hand, from a theoretical point of view, 
this work represents a twofold extension of the math- 
ematical treatment originally introduced by Huber. [2(J 
In the first place, instead to study the asymptotic fluores- 
cence decay in random Frenkel lattices we consider the 
case of general aperiodic systems. In the second place, 
meanwhile Huber evaluated the dipole moment matrix 
elements using unperturbed eigenvectors and perturbed 
eigenvalues, we show that this approach can only explain 
the occurrence of long-period oscillations. To account for 
sho rt-pe riod oscillations we must include the second term 
of (10a), that is, perturbed eigenvectors. This second 
term is responsible for the coupling between the lowest- 
lying mode and the high-lying modes through the aperi- 
odic arrangement of traps, causing the short-period os- 
cillations. 



VI. CONCLUSIONS 

In this paper we have carried out a detailed analysis of 
the fluorescence decay in self-similar Frenkel lattices with 
interstitial traps arranged according to the Fibonacci and 
Thue-Morse sequences. To this end, we have solved nu- 
merically the equation of motion of Frenkel excitons in 
the time domain. Unlike what it was previously found in 
random lattices, we have observed that the fluorescence 
decay after <5-pulse excitation, in which all centers are 
initially excited, behaves exponentially on average in the 
case of aperiodic systems. We have demonstrated that 
such a time decay is quite well described within a ho- 
mogeneous medium approximation, where the aperiodic 
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distribution of traps is replaced by a uniform distribution 
one with a renormalized trapping rate. This exponential 
trend is similar to that found in the case of incoherent 
transport through aperiodic lattices and, in spite of its 
dependence on the concentration of traps, it cannot be 
used to fully characterize each particular arrangement of 
trapping centers. Instead, quantum effects reveal them- 
selves through small-amplitude oscillations superimposed 
to the average exponential decay. Interestingly, we have 
shown that the oscillatory pattern is characteristic of 
each particular lattice and, consequently, it can be con- 
sidered as a fingerprint that allows us to elucidate the 
kind of underlying ordering present in the distribution of 
traps. 

We have carried out a perturbation approach to ac- 
count for the small-amplitude oscillations, finding an ex- 
cellent agreement with numerical data. The perturba- 
tive analysis clearly indicates that the oscillatory pat- 
tern of the fluorescence decay is directly related to the 
Fourier transform of the distribution of traps and, as a 
consequence, it is characteristic of the ordering present 
in the considered system. In this sense, we realize that 
the necessary condition for an arbitrary system to ex- 
hibit Fourier peaks is the occurrence of long-range or- 
der in its structure. This condition plays a major role 
in our theoretical analysis. Since each well-defined fre- 
quency contributing to the function p(t) can be related to 
a distinctive peak in the Fourier transform of the lattice, 
the Fourier analysis of the time-dependent signal gives 
us direct information concerning the spatial distribution 
of traps. The physical origin of the fluorescence oscil- 
lations can be traced back to the coupling between the 
low-lying exciton mode and certain particular high-lying 
modes through the lattice topology. This means that ex- 
citons act as a probe of the spatial distribution of trap- 
ping centers. Hence the detailed analysis of time-decay of 
the fluorescence, which uses excitons as a physical probe 
of the long-range order, represents a powerful tool to in- 
vestigate the topology of general aperiodic systems from 
experiments, opening in this way a novel experimental 
technique which may be tentatively referred to as laser 
crystallography. 
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